DISTANCE DETERMINATION OF VARIABLE 
GALACTIC SOURCES 



J. Hu\ S. N. ZHANGi'2,3,4 rj. p ^^1,4 

^ Physics Department and Center for Astrophysics, Tsinghua University, Beijing 100084, 
P.R. China 

Physics Department, University of Alabama in Huntsville, Huntsville, AL ,35899, USA 
Space Science Laboratory, NASA Marshall Space Flight Center, SD50, Huntsville, AL 
35812, USA 

^Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, China 



Abstract 



We have developed a timing analysis method to determine the dis- 
tances of variable galactic X-ray sources based on the method advanced 
by Triimper and Schonfelder in 1973. The light-curve of the halo pro- 
duced by the scattering of X-rays off the interstellar dust is delayed 
and smeared by the dust grains. This method utilizes the differences 
between the power density spectra of the point source and the halo. 
We present the details of this method and our first applications of this 
method to the Chandra data of X-ray binary Cyg X-3. 
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1. Introduction 

When the high energy photons emitted by the X-ray source arrive 
to the earth, they have been scattered by the interstehar dust grains. 
Therefore many galactic X-ray sources are observed to have faint and 
diffuse halos. This effect and its possible application was first discussed 
by Overbeck (1965), and the X-ray halo was first detected with Einstein 
Observatory (Rolf 1983, Catura 1983). The investigation of the shape 
and strength of the haloes was done by several authors (e.g. Mauche & 
Gorenstein 1986, Predehl & Schmitt 1995) based on the data of Einstein 
and ROSAT. 
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As the photons in the halo travel through longer paths than the pho- 
tons in the point source, any intensity variation in the the source in- 
duces a delayed signal in the halo. It has been suggested by Triimper 
and Schonfelder (1973, hereafter TS1973) that this effect can be uesd to 
determine the distance of variable galactic X-ray sources. Although the 
method is very simple, the first successful application was realized by 
Predehl (2000) for a X-ray binary Cyg X-3 almost 30 years later. They 
obtained the distance of ^V"^ to the source. The difficulty of the 
method is that it is hard to distinguish the halo from the image induced 
by the point spread function of the instrument until the Chandra X-ray 
Observatory provides us the high resolution image of X-ray sources. 

The shape and strength of the halo depend on the distance to the 
X-ray source, the geometric distribution and the physical properties of 
the interstellar dust grains (Hayakawa (1970), Mathis & Lee (1991), and 
Predehl k Klose (1996)). 

2. Methodology 

The differential scattering cross section of X-ray photons in interstellar 
dust is (Predehl, 2000) 

^ = |5($)|2 oc exp(-0.46£;2a2$2) 

where a is the size of the dust grains (in ^um) , E is the energy of X-ray 
(in keV), and $ is the scattering angle (in arc min). 




Figure 1. Scattering geometry. indicates the observer, d indicates the source, xd 
is where the scattering occurs. 



Following TS1973, we obtain the time delay of a scattered photon 
compared with an unscattered photon is (See Fig. 1) 

t = —- = 1.15d- , (2) 

2cl-x 1 -x' ^ ' 

where d is given in kpc, G in arcsec. 
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Given the dust density distribution along the hne of sight p{x), the 
time delay distribution function H{t, @) can be calculated. As pointed 
by Predehl & Schmitt (1995), we can assume p{x) is a constant. Thus 

exp(-i^[(l + -^)f]2) 

For a viable point source with light-curve I{t), the normalized light- 
curve in the halo is 

r/ , I(T)H(t - T, @)dT 

where tm(0) is the possible maximal time delay. 

In TS1973, I{t) and B{t, B) are compared directly to determine dis- 
tance. 

Suppose the light-curve I{t) of the point source and the light-curve 
B(t,oj) of the halo at uj radius are stationary random signals, and 
h{u),Q) are their Fourier transformation respectively. It is easy to show 
that 

h{uj,Q) ^i{uj)h{uj,Q) (5) 

where 

h{LO, 0) = y H{t, e)e-*'^*di ^ J H{t, e)dt (6) 

is the normalized integral transformation of the signal arrival time dis- 
tribution function H{t,0). 

Denote the power density spectrum of the point source and the halo 
as Pb{io, Q) and Pi{co), we get 

Pb{u;,e) = Pi{Lo)h\u;,e) (7) 

Given the power density spectra of the point source and the halo at 
different radius, h{uj,@) and H{uj,Q) are easy to know, and d can thus 
be determined by fitting equation (3). 

3. Application to Chandra data on Cyg X-3 

Now we apply this method to the Chandra data of X-ray binary Cyg 
X-3. Cyg X-3 is an eclipsing X-ray binary with an orbital period of 
4.8 hours (Brinkman et al. 1972). Cyg X-3 was observed with the 
High Energy Transmission Grating Spectrometer (HETGS) on board of 
Chandra on October 20, 1999, with a total time of 12.3 ksec, starting 
at 01:11:38 UT. The observation covered the binary phases from -0.3 to 
-1-0.4. We only use the zeroth order data in this study. 
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Figure 2. The theoretical value of the Function h'^{LO,€> = A")/h?{uj,Q = 2") at 
5.0-7.0 energy band. The three curves take d = 5, 9, 13 kpc respectively. 



Because the image of the point source is severely damaged by the 
pileup effect, we can only campare the power density spectra at different 
and radius. The theoretical difference of the power density spectrum of 
the halo at 2 arcscc and 4 arcscc in 5.0-7.0 keV energy band are shown 
in Fig. 2. Then we calculate h'^{LJ,Q = [4" , 6"]) / h'^ {to , @ = [2", 4"]) in 
5.0-7.0 keV energy band, the result shows in Fig. 3. Because the error 
bar of the data in Fig. 3 is very big, it is hard to fit the data with 
theoretical curves. For simplicity, we fit the data with a exponential 
curve. The turning point of the fitting curve lies at about 0.001 Hz, and 
the corresponding value of d is 9 kpc. According to Fig. 3, it seems the 
real turning point should small than 0.002 Hz, so we estimate d is at 
least 5 kpc. 
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Figure 3. Observed value of the function h'^{u;, 6 = [4", 6"])/h^{Lu, & = [2", 4"]) in 
5.0-7.0 keV energy band. 
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4. Discussion 

Our result is consistent with previous estimate on the distance of Cyg 
X-3. Predehl et al. (2000) used TS1973 method and obtained 9ti kpc, 
Dickey (1983) found a lower limit of 9.2 kpc 1 using 21 cm wavelength 
absorption data, Predehl &; Schmitt (1995) derived 8 kpc as distance 
through the galactic dust layer from their comparison of X-ray scattering 
and absorption. 

Our method and TS1973 method both determine the geometric dis- 
tance instead of physical distance. The physical distance depend on the 
estimation of many physical parameters, i.e., absorption optical depth. 
The uncertainties of the geometric distance determination come mainly 
from the dust density distribution along the line of sight. In principle, it 
can be solved by carefully fitting the shape and strength of the halo. In 
TS1973 method the X-ray source must have relatively obvious variations 
in the time interval comparable to the delay time (see Fig. 4. in Predehl 
et al. 2000). But in our method this requirement is not needed. 

Although the Chandra Observatory has high resolution, its effective 
area is relatively small, therefore the count rate of the Cyg X-3 halo is 
not high enough for good statistical study with our method. The future 
X-ray satellites such as Constellation- X will provide us with high count 
rate data, and our method can be applied to many more galactic variable 
X-ray sources. 
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